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The nonlinear interaction, due to quantum electrodynamical (QED) effects, between an electro¬ 
magnetic pulse and a radiation background is investigated, by combining the methods of radiation 
hydrodynamics with the QED theory for photon-photon scattering. Eor the case of a single coherent 
electromagnetic pulse we obtain a Zakharov-like system, where the radiation pressure of the pulse 
acts as a driver of acoustic waves in the photon gas. Eor a sufficiently intense pulse and/or back¬ 
ground energy density there is focusing and subsequent collapse of the pulse. The relevance of our 
results for various astrophysical applications are discussed. 
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The theory for electromagnetism in vacuum can be 
roughly divided into that for coherent and that for in¬ 
coherent systems. In the former case Maxwells equations 
are the explicit starting point, and that theory leads 
to most of the wellknown applications. In the latter 
case, the ensemble average of the electromagnetic field 
is zero everywhere, and the system can then essentially 
be treated as a gas of ultra-relativistic particles 0 . Most 
of the theory for radiation gases cannot be considered as 
pure “electromagnetism in vacuum”, since the interac¬ 
tion of photons with other particles is one of the main 
features of the theory. In the present Letter, we will con¬ 
sider the self-interaction between photons, due to quan¬ 
tum electrodynamical (QED) effects. One of the many 
interesting aspects of such an analysis is the possibility 
of photon-photon scattering, due to the interaction of 
light quanta with virtual electron-positron pairs 0- By 
integrating out the virtual pairs, one obtains an effective 
field theory for the photon-photon interaction, in terms 
of the electromagnetic field variables (E/u/). The lowest 
order correction to Maxwell’s vacuum equations is then 
conveniently expressed by means of the Heisenberg-Euler 
Lagrangian [see Eq. below]. In that way many pa¬ 
pers on QED photon^hoton scattering have been writ¬ 
ten (e^.j_0 0 0 IE 0 0 0 0 and references therein). 
Refs. 00 I 3 concern techniques for laboratory detection 
of QED effects, involving second harmonic generation 0 , 
self-focusing [ffj, and nonlinear wave mixing in cavities 
0, respectively. Other aspects of general theoretical in¬ 
terest that have been dealt with are, e.g., effects of static 
magnetic background helds on higher harmonic genera¬ 
tion, wave propagation velocities in the QED vacuum and 
one-loop corrections to the Heisenberg-Euler Lagrangian 
0 , two-dimensional collapsing scenarios 0 , and the re¬ 
fractive properties of the QED vacuum |9j. However, 
the above studies involve strictly coherent helds, or one- 
photon phenomena. For instance, the non-trivial propa¬ 
gation of photons in strong background electromagnetic 
helds, due to effects of nonlinear electrodynamics, has 
been considered in a number of papers (see, e.g.. Refs. 


0 0 and references therein). The main focus in these 
papers was on the interesting effects of photon s plit ting 
and birefringence in vacuum. However, Thoma |iri in¬ 
vestigated the interaction of photons with a photon gas, 
using the real-time formalism, and calculated the corre¬ 
sponding change in the speed of light due to the Cosmic 
Microwave Background (CMB). 

In the present letter we are going to combine methods 
from radiation hydrodynamics 0 with QED theory for 
photon-photon scattering, in order to give a framework 
for the interaction of coherent electromagnetic helds with 
a radiation gas background. For the case of a single co¬ 
herent electromagnetic pulse we obtain a Zakharov-like 
system EH , where the radiation pressure of the pulse acts 
as a driver of acoustic waves in the photon gas. Similarly 
to ordinary acoustic waves these waves are longitudinal 
and characterised by variations in pressure, and thereby 
in energy density. The index of refraction depends on 
the photon gas energy density, and thus the excitation of 
the acoustic waves leads to a back-reaction on the pulse. 
For a sufficiently intense pulse and/or background energy 
density we may have focusing and subsequent collapse of 
the pulse. Applications to coherent pulses propagating 
in the present as well as in the early cosmic radiation 
background are discussed, together with other astrophys¬ 
ical phenomena. While there are possible explanations 
for the small CMB structures recently detected , we 
think that our mechanism below for CMB structure for¬ 
mation involving photon-photon interaction should be of 
interest. 

The nonlinear self-interaction of photons can be for¬ 
mulated in terms of the Heisenberg-Euler Lagrangian 0 

L = eoF + Kel -k TG^] , (1) 

where F = (F^ — cPB^)j2 and G = cE ■ B. Here 
K = 2a^Ti^/ k, 1.63 x 10“^°ms^/kg, a is the 
fine-structure constant, h is the Planck constant, rrie the 
electron mass, and c the velocity of light in vacuum (for 
higher order corrections, see, e.g., expression (26) of Ref. 
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0). The Lagrangian iPl is valid as long as there is no 
pair creation and the field strength is smaller than the 
critical field, i.e., 

w ^ /fi and \E\ <C i?crit = me(? jeXc (2) 


where r denotes the group velocity of the photon, k the 
force on a photon, and the dot denotes time derivative. 

The equation for the collective interaction of photons 
can then be formulated as M 


respectively. Here e is the elementary charge, Ac is the 
Compton wave lengtln and Edit — 10^® V/m. 

According to Ref. Q , the dispersion relation for a low 
energy photon in a background electromagnetic field can 
be derived from the Lagrangian o, with the result (see 
also Ref. 0 and references therein) 

Loik,E,B)=c\k\{l-XX\Q\^). (3) 


(9) 

where the distribution function has been normalised such 
that the number density is n{r,t) = f N{k,r,t)dk. 

For a general function f{k,r,t), the moment equation 
is 


where 
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and A = A±, where A+ = 14k and A_ = 8k for the two 
different polarisation states of the photon. Furthermore, 
k = k/k. The approximation AjQp <c; 1 has been used. 
The background electric and magnetic fields are denoted 
by E and B, respectively. 

We will below study two scenarios: (i) A plane wave 
pulse propagating on a background consisting of a radi¬ 
ation gas in equilibrium, and (ii) a radiation gas affected 
by an electromagnetic (EM) pulse propagating through 
the gas. For case (i), the relations (fep • E)'^ = ^E^, 
{kp ■ B)^ = and E x B — 0, hold, where kp is the 
pulse wave vector. Hence, from 0 we obtain 

IQgasi' = 1^, (5) 

where S’ = eo{E^ + c^i?^)/2 is the energy density of the 
radiation gas. For case (ii), we argue in a similar manner. 
To lowest order, the directions k of the photons in the 
gas are approximately random, and the EM pulse is a 
superposition of uni-directional plane waves, such that 
E = EpB, and B — Epkp x ejc, where e is the unit 
electric vector. Then 0 yields 

IQpuisel" = leolEpl'^. (6) 

With the relation 0 , and using standard methods for 
slowly varying envelopes m , we then derive the dynam¬ 
ical equation for a pulse on a photon gas background 

i -|- ckp ■ Ep + '^XckpSEp = 0, (7) 


where — {kp ■ V)^, noting that the self¬ 

interaction of the pulse vanishes. 

For a dispersion relation uj = ck/R{r, t), where R is the 
refractive index, we have the Hamiltonian ray equations 


r 
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duj c f 
-Vo; = |vi?. 


(8a) 

(8b) 


where (/) = n{r,t) fNdk. Choosing / = huj, we 
obtain the energy conservation equation 


dS 

'm 


V • {S’u q) 


R~m’ 


(lla) 


where S{r,t) = nfhcu) is the energy density, and q{r,t) = 
n {hujw) the energy flux. Here we have made the split 
r = u w, where {w) = 0. Thus w represents the 
random velocity of the photons. With / = tik, we obtain 
the momentum conservation equation 
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(11b) 


where H = n {fik) is the momentum density, and 
P = n {w ® {kk)) is the pressure tensor. It follows im¬ 
mediately from the definition of the pressure tensor that 
the trace satisfies Tr P = n {Tikw ■ k) = nR {Tllow ■ k)/c. 
For an observer comoving with the fluid, i.e., a system in 
which (tt)o = 0 (the 0 denoting the comoving system), 
Eq. 18all shows that {w ■ fc)o = {R)o/c, so that the trace 
of the pressure tensor in the comoving system becomes 
(TrP)o = (^)o- For an isotropic distribution function, 
the pressure can be written as P = Tr P/3, satisfying the 
equation of state P = Sf3. 

The system of equations llllall and (HTbll needs closure, 
and this can be obtained by choosing an equation of state. 
Furthermore, we still have some gauge freedom since we 
have not specified the frame. Two choices prevail in the 
literature, (i) the energy frame, in which q = 0, and (ii) 
the particle frame, where u = 0. We will here adopt 
a particle frame, and furthermore choose the equation 
of state to be P^- = SijSl3. The relation between the 
heat flow and the momentum density is straightforward 
to derive. The result is q = c^II/R^. From the disper¬ 
sion relation 0 and Eq. we have i? ~ 1 + |Aeo|PpP, 
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from which, through Eqs. CH), we obtain 


_+c V-JT - -Aso 




dt 


+c2v • (iT|Ep|2) 


,(12a) 


and 


+ 3^^ “ gAeo<^V|£'pp. 


(12b) 


We now introduce ^ = <^o + where <fo 5S) is 
the unperturbed energy density in the absence of the elec¬ 
tromagnetic pulse, noting that for the unperturbed state 
11 = 0. Taking the divergence of 112bll and eliminating 
n using 112 all . we obtain 
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_V <5^ _ ^ 


\Ep\^ + 


\Ep\^ 

dt^ 


(13a) 

describing acoustic waves in a photon gas driven by a 
coherent electromagnetic pulse. We note that the set 
of equations 113all and 0 resembles the wellknown Za¬ 
kharov system, describin g co upled Langmuir and ion- 
acoustic waves in plasmas |l lj . Transforming to a system 
moving with the pulse r = t, ^ = x — ct (letting the pulse 
wave-vector be kp = kpx), we have 




— — — Aeo^o 
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d^\Epf 

oe 


(13b) 

dropping the slow derivatives proportional to dfdr. For 
a general pulse geometry Eq. OB requires extensive 
analysis, but for specific forms of the pulse we can sig¬ 
nificantly simplify the expression for 6S'. For a broad 
electromagnetic pulse where V_l <C d/d^ we thus inte¬ 
grate H13a|l to obtain SS" = —2XeoS’o\Ep\‘^, whereas for 
a long needle shaped pulse with d/d^ <C we instead 
have 6S' = 2XeoS'o\Ep\‘^. The unperturbed photon gas 
energy density S’o will give a frequency shift for the solu¬ 
tion of Eq. 0 , but this shift can be transformed away. 
The resulting equation is 


riW r A 

T -X^ckpe^<g^\Ep\^Ep = 0, (14) 


where the minus (plus) refers to a broad (needle shaped) 
pulse. For a broad pulse, wave collapse does not oc¬ 
cur, but for a needle shaped pulse, collapse will occur 
when the nonlinear term dominates over the diffraction. 
Comparing the terms in Eq. d, we can obtain a rough 
estimate of the collapse criterion, namely 


(15) 


which is consistent with Ref. M Here Vp denotes the 
pulse width, S'o = eol^^gasP, and |Echar| = (eo't)”^^^ « 
2.6 X 10^°V/m. Although the collapse of electromag¬ 
netic pulses is complicated la, the qualitative fea¬ 
tures can be described as a divergence in the pulse en¬ 
ergy density, while its width decreases to zero, in a 
finite time. A rough estimate of the collapse prop¬ 
erties can be obtained assumi^ cylindrical symme¬ 
try of the pulse. Following [l^, the trial function 
Et{t, r) = A{t) sech[r/a(T)] exp[j5(r)r^] is used together 
with Rayleigh-Ritz optimisation, in order to reduce the 
problem to a differential equation for the width a{T). The 
field strength then satisfies |A(t)|/|A(0)| = a(0)/a(T), 
and the width behaves as (a(T)/a(0))^ — 1 oc (1 — 
a(0)^|A(0)p//c)T^ where Ic — 0.5 x (fc^A^eo'^o)”^ Thus, 
if a(0)^|A(0)p > Ic the pulse will collapse to zero width 
in a finite time. Of course, when the pulse intensity in¬ 
creases, higher order effects will become important, possi¬ 
bly halting the collapse. Furthermore, as the pulse width 
decreases, derivative corrections to o become important 
113 , and thus change the collapse scenario. 

The most intense pulses in our universe are the gamma- 
ray bursts Powers of the order 10“^^ W in the 

gamma range (kp ~ 10^^ m“^) means that we may have 
\Ep\kprp ^ 10^® V/m [^. In todays universe the energy 
density of the CMB is 2.6 x 10® eV/m®, which correspond 
to |Agas| ^ 7 X 10“^V/m. Thus the collapse condition 
is (far from surprising) not fulfilled for gamma-ray 
bursts propagating in the present microwave background. 
However, we note that if there were mechanisms for gen¬ 
erating equally intense gamma pulses in the early uni¬ 
verse at an age t ^ 6 x 10 ® yr, such pulses should indeed 
have collapsed, as for those times the energy density of 
the CMB obeyed |Fgas| ^ 2 x 10"^ V/m implying that itTB 
is fulfilled. 

Furthermore, we note that collapse may occur for much 
less intense pulses than gamma-ray bursts, in case the 
propagation takes place in a photon gas which is more 
energetic than the CMB. Highly energetic photon gases 
should exist today in the vicinity of pulsars and mag- 
netars. For pulsars, the low-frequency dipole radiation 
is not able to leave the system directly, as it is first re¬ 
flected by the surrounding plasma. The momentum of 
the low-frequency photons in the vicinity of pulsars are 
not necessarily thermally distributed, but still we should 
be able to think of the pulsar environment as a highly 
energetic photon gas. For magnetars (i.e., pulsars with 
surface magnetic field strength ~ 10^° T |lJl), the high 
fields create a surface tension leading to star quakes, and 
subsequent generation of low-frequency photons. This 
environment could be in an electromagnetically turbu¬ 
lent state [i^. Thus electromagnetic pulses of a com¬ 
paratively moderate intensity propagating in such envi¬ 
ronments may collapse, due to the high energy densities 
of the photon gases. Rough estimates give energy densi¬ 
ties of such photon gases that can be 20 — 22 orders of 


kyp\Ep\^\E,^\^ > \Ec^y 
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magnitude larger than that of todays CMB. 

While the conditions for collapse may be fulfilled only 
during extraordinary circumstances, there is still a possi¬ 
bility that energetic pulses, although not intense enough 
for collapse, can leave a certain inprint in the CMB. Since 
we are mainly interested in effects that persist after the 
EM-pulse has passed a given area, we need to abandon 
the approximation of a broad or needle shaped pulse and 
instead investigate solutions to 11 .Ihll . If the radiation gas 
energy density is well below that for collapse, the back- 
reaction on the EM-pulse is small, and we do not have to 
solve the system 0 and (Ebl) self-consistently, but can 
take the EM-pulse as given. Furthermore, investigating 
the generation of acoustic waves during times short com¬ 
pared to the diffraction time, we can consider the pulse to 
be static in the comoving frame. To simplify the analysis 
by making the r and ^-dependence separable, we assume 
that the source term obeys V5_ \Ep\^ + \Ep'f' = 0, 
where rg is the characteristic radius of localisation. For 
a pulse with no angular dependence, this means that we 
can write \Ep'f' = Jo{r/ro)S{^), where Jq is the 0*^*^ order 
Bessel function, and we leave S{^) unspecified. Letting 
the generated acoustic waves have the same radial de¬ 
pendence as the pulse, i.e. 6S‘{r,^) = Jo{r/ro)£’^{^) we 
then obtain from 11 ,'Ib|l 


de 




Xsof^o 



2 


d^S ' 
d^\ ■ 


(16) 


With the boundary condition of no acoustic waves before 
the pulse has passed, the solution of m is 


<§’5(0 = SiO sin[(^ - a/V2ro] d^'. (17) 

^0 J oo 

Apparently, after the pulse passage the field can be writ- 
ten as sin[(^/-\/2ro) -I- <5] , where the ampli¬ 

tude and the phase angle 5 depend on the detailed 
form of the pulse profile 5(0) (see Ref. for analyt¬ 
ical expressions). For a pulse length slightly less than 
ro we have ^50 ^ Aeo<^o5max ~ Aeo<§’o where 

lEpl^ax central value of \Ep'^. Currently, mea¬ 

surements of the CMB can detect relative temperature 
anisotropies of the order of 10“® Thus, we see that 
l^'pl 0 10 ^'^V/m must be fulfilled for possible detection 
of the resulting background anisotropy. Our estimate for 
the acoustic wave amplitude applies only relatively close 
after the pulse passage, before the acoustic wave begins 
to spread due to diffraction (as associated with the slow 
T-dependence omitted in Eq. iTT^ I. It is clear from the 
above that intense pulses leave an inprint in the CMB 
after the pulse has passed a given area, and the results 
shows that anisotropies in the high-frequency electromag¬ 
netic spectrum partly transfer to the lowfrequency back¬ 
ground. However, our calculations apparently do not in¬ 
clude all effects that will influence the earth-based mea¬ 


surements. Thus much work is still needed in order to 
determine whether the effects induced in the CMB by 
astrophysical sources can be seen in the detailed experi¬ 
ments that are currently made |22l| . 
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